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In this paper we present a model of a symmetric Brownian motor (SBM) which changes the sign 
of its velocity when the temperature gradient is inverted. The velocity, external work and efficiency 
are studied as a function of the temperatures of the baths and other relevant parameters. The motor 
shows a current reversal when another parameter (a phase shift) is varied. Analytical predictions 
and results from numerical simulations are performed and agree very well. Generic properties of 
this type of motors are discussed. 



I. INTRODUCTION 

We all do know that it is possible to extract some 
amount of mechanical work from a thermal bath at a 
temperature Ti provided we have another bath at a lower 
temperature T% < Ti. Thermal engines are the devices 
that perform this task. All this is well known from el- 
ementary textbooks on thermodynamics. We also know 
from statistical mechanics that any object in a thermal 
bath exhibits random energy fluctuations of the order 
KbT. These fluctuations are relatively very small for 
macroscopic objects but of very important relevance for 
nanometric objects such as biological motors: kinesins, 
dyneins, etc. We are also familiar with windmills, 
which are able to extract useful work from random winds 
by a proper adaptation to the wind direction. We can ask 
ourselves if it is possible to rectify thermal fluctuations 
by some appropriate mechanical devices. 

The engines which aim to get useful work by rectifying 
thermal fluctuations are called Brownian motors (BM). 
In fact, the paradigm of such speculations is Feynman's 
famous ratchet anS pawl machine 0. During the last 
years a lot of effort has been invested to study what 
has been called as the ratchet effect. This is a mech- 
anism which consists in breaking the spatial and tem- 
poral inversion symmetry of the system so that directed 
transport emerges, often enhanced by the thermal fluc- 
tuations. The ratchet mechanism can be implemented 
in different ways. Here we will make a model through 
an equation for a dynamical classical variable (position 
or angle) moving in a periodic and asymmetric poten- 
tial (a ratchet potential) coupled with another degree of 
freedom which will break thermal equilibrium. 

There is an enormous variety of ratchets in the liter- 
ature: pulsating ratchets (on-off, fluctuating potential, 
traveling potential), tilting ratchets (fluctuating force, 
rocking ratchets), Seebeck ratchets, Feynman ratch- 
ets, temperature ratchets, frictional ratchets, quantum 
ratchets, collective ratchets, mechanico-chemical ratch- 
ets, ionic bombs and pumps, etc. |3(. 

In this work we will focus on a new Brownian mo- 
tor inspired by the so called Feynman's ratchet (Fig^). 
Feynman stated in his famous lectures Q that a very 
particular device (the ratchet and pawl machine) could 




Figure 1: Feynman's original ratchet and pawl machine. 



Figure 2: Detail of the stopping mechanism based on the 
interaction between the ratchet and the pawl. 



have the efficiency of a Carnot engine when operating 
reversibly. He idealized a gadget (see FigQJ made up 
of two boxes with some gas at different temperatures. 
The hotter box contains an axle with vanes in it. The 
bombardments of gas molecules on the vane make the 
axle rotate with random symmetric fluctuations. At the 
other end side of the axle there is a second box with a 
toothed wheel which in principle can turn only one way. 
The pawl (the stopping mechanism) is under the influ- 
ence of another temperature (Fig|2J) • At first glance one 
could think that it seems quite likely that the wheel will 
spin round one way and lift a weight even when both 
gases are at the same temperature and thus violating the 
Second Law. However, a closer look at the pawl reveals 
that it bounces and so the wheel will rotate randomly 
in any direction, doing a lot of jiggling and with no net 
turning. Thus, the machine cannot extract work from 
two baths at the same temperature. 

When the temperature of the vanes is higher than the 
temperature of the wheel, Feynman concludes that some 
work is performed with Carnot 's efficiency when the ma- 
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chine is lifting the weight very slowly. This is indeed a 
very optimistic result which has been revised, many years 
later, in Refs. 0,01. Indeed, there are some overlooked 
aspects in Feynman's argument which can be summa- 
rized as follows: since the engine seems to be simultane- 
ously in contact with two baths at different temperatures 
(through the rigid axle), it cannot work in a reversible 
way and Carnot's efficiency will never be achieved. We 
will come back to this point along this paper. 

A possible mathematical model of the device of FigQ] 
and FigElin terms of overdamped Langevin equations is, 
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where £j (t) mimics thermal fluctuations and it is assumed 
to be a white noise, Gaussian distributed with zero mean, 
satisfying the fluctuation-dissipation theorem, 

< >= 2fc B T i A i $(t - t')8 lv (3) 

The total potential energy is modeled by, 

V{6,h)=T6+\kh 2 + , (4) 

e 'o — 1 

where, Hr{&) represents the periodic but asymmetric pro- 
file of the ratchet. An explicit expression will be given in 
the next section. 

The potential V has three terms, t is the torque which 
gives useful work. The second term accounts for the po- 
tential energy stored in the spring (the pawl) with con- 
stant k, which pushes down the ratchet. Finally, the last 
term is a very repulsive potential at short distances and 
nearly zero at long ones. It is used to avoid that the pawl 
crosses through the real physical surface of the ratchet. 
At the same time, it couples both degrees of freedom 
and h. Other potentials can account for the same physics 
in this mathematical scheme 0, |(| . 

Numerical simulations show a non-zero mean velocity 
of the ratchet device for very high temperature differ- 
ences. The main conclusion is that, although the motor 
can perform useful work, it has an extremely small effi- 
ciency, being very far away from Carnot's efficiency es- 
timated by Feynman, as pointed out and shown in Refs. 
m IS 111 • This motor presents a high thermal conduc- 
tivity and as a consequence has a very low efficiency. 

The structure of this paper is the following. In sec- 
tion II we present the model for a symmetric Brownian 
motor and the numerical results obtained by computer 
simulations of the equations of such device. Section III 
is devoted to the analytical approach to this model. We 
end in Section IV with some comments and conclusions. 
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Figure 3: Model for the symmetric Brownian motor (SBM). 



II. A SYMMETRIC BROWNIAN MOTOR 
(SBM) 

Feynman's ratchet, and similar models, do not fulfill 
the following inversion property: T 2 «-> T\ — ► v —v. 
Thus our initial motivation to propose this new model 
is to have a motor such that, when switching the tem- 
peratures of the two baths, this is, when reversing the 
temperature gradient, the mean velocity absolute value 
does not change while the sign of the velocity does. Ob- 
viously, this would imply some geometrical symmetries 
in the engine. One example of a microscopic model for 
a Brownian motor that suffers this change in the sign of 
the current can be found in [8j. There is also a model 
for a diode rectifier that is analytically solvable and 
that has the same main features as the SBM, namely, it 
is symmetric and it shows " reversibility" . 

Our "symmetric" motor is shown in Fig [31 According 
to this scheme, the stochastic differential equations that 
describe the dynamic evolution of the system in the over- 
damped regime are, 

M^- = -de 2 V + b(t). (5) 

The noise terms are again white thermal noises, Gaus- 
sian distributed and with zero mean, satisfying the 
fluctuation-dissipation theorem yj. 

The potential V, contains an external torque r, two 
ratchet-shape potentials and an harmonic interaction 
that couples both degrees of freedom, 

v(e 1 ,e 2 ) = t0i + Vr(0i) + v R (-9 2 ) + h(e 1 + <t>- e 2 )\ 

(6) 

where the phase shift <f> is one of the most important 
parameters. 

In FigQ]we see the form of the ratchet potential Vr{6) 
used in the numerical simulations. Its analytical expres- 
sion is, 



V R (6) 
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(7) 



V(0) = sin(d0) + 0.275 sin(2d(9) + 0.0533 sin(3d0). (8) 
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Figure 4: Shape of the ratchet potential Vr for d = 4. 



Vo controls the height of the potential, d is the number of 
teeth per cycle 27r, and the asymmetry of the potential 
is controlled by changing the numerical coefficients that 
multiply the sinus functions expansion of the potential. 

Notice that the ratchet potential that one variable (8i) 
sees is the specular image of the one that the other vari- 
able (02 ) feels, as in the diode rectifier of Ref. 0- It is 
in this sense that we call our model symmetric. 

To study the relevance of each parameter we proceed 
to express our system in terms of dimcnsionless ones. In- 
troducing the dimensionless time s as t = ^s, the set 
of equations JSJ) becomes more compact. The dimension- 



less parameters are now: T\ — 



_ K B T 2 



U — k 



Vb ' - 1 * Vo 

A = j^, t = <pr and, of course, <f>, which is in radians. 
Notice that Vq controls the energy scale, which is mea- 
sured in KbT\ units. We also see, for instance, that only 
the fraction of the friction coefficients is relevant. 
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Figure 5: Mean velocity versus T2/T1 (r = 0). 
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Figure 6: Mean velocity of the motor at the optimal regime as 
a function of the external torque r. The stall torque (where 
v — 0) is around r = 1.2. 



Numerical results of SBM 



Preliminary numerical results indicate that the motor 
inverts its velocity when the temperature gradient is also 
inverted (see Fig[5J). Some important thermodynamic re- 
quirements are also fulfilled: there is no net motion in 
the limits Ti — ► T\ (a unique bath), k — > (no coupling 
between baths) and k — » 00 (only one effective tempera- 
ture) . 

To perform the simulations we take A = 1 for simplic- 
ity and KbT\ = 1. As in previous models for motors 
we have explored the parameter domain to get the most 
effective values (larger velocities). This situation corre- 
sponds to the parameter values Vq = 2.5, k = 100 and 
d = 16. In FigOand Fig^Jsome numerical results of the 
mean velocity v = (9\) are presented for = 0.4. We see 
in Fig [3] that there is a maximum around TijT\ = 2.5 
and a minimum at T2/T1 = 0.5. This implies that the 
velocity is bounded and larger temperature gradients do 
not imply larger velocities. For very large T2 the motor 
does not see the ratchet potential and cannot take advan- 
tage of the broken inversion symmetry. Then the motor 



has a zero mean velocity. The linear dependence of the 
velocity around T%jT\ ~ 1 is is clear signature that the 
motor has the inversion property. We will extend this re- 
sult further. In FigElwe see the expected negative linear 
dependence of the velocity versus torque until the stall 
force where the velocity is zero. 

Moreover the most striking result is the velocity inver- 
sion as a function of the phase difference parameter <fi. 
We show in FigEI the mean velocity of the motor as a 
function of the phase shift 4> in the cases of T2 = 2Ti and 
when T\ = 2T2. One must notice that both cases are 
" symmetric" . This result allows to extract some conclu- 
sions: the existence of a current inversion as a function of 
4>, the great amplification of the speed for the particular 
value of <p = 0.4 and, finally, that the motor is perfectly 
symmetric when reversing the temperature gradient. 

These results deserve theoretical explanations that will 
be presented in section III. 
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Figure 7: Current inversion profile v(<f>) when T2 = 2T\ (cir- 
cles) and when Ti = 2T2 (triangles). The angle goes from 
to 0.52 because we are considering d = 12 and therefore the 
periodicity is || = 0.52. The dashed lines are a guide to the 
eyes. 



B. Energetics of the SBM: the efficiency 

It is not enough to know whether a motor will move, 
nor even how fast it will run. One would also like to 
know how efficiently such machine can operate. There 
are many different approaches in the literature to study 
energetics and, in particular, efficiencies [EEIE!- We 
will follow Sekimoto's characterization of the energetics 
of thermal ratchet motors @ . This scenario is the most 
suitable scheme both because it is intuitive and also from 
its potential applicability. The main idea is to find how 
much energy, received from an external source, the mo- 
tor will be able to employ to produce useful mechanical 
work, i.e. to lift a load. In order to compute this, we 
must understand where the energy comes from, where it 
can go and what ingredients account for that. Sekimoto 
deals with energies and not with power. When one tries 
to study the problem in units of power, velocities ap- 
pear explicitely in the computations and, since they are 
instantaneously not continuous because of the stochas- 
ticity of the system, one finds that some of the terms 
are extremely bad behaved and it is impossible to get 
a reasonable numerical convergence in the simulations. 
Moreover, we have found it easier to directly compute 
the stochastic integrals and afterwards perform the aver- 
age. 

Let's analyze the energetics of such kind of systems. 
We would expect there's some energy R that is released 
from the bath at T 2 (the one at a higher temperature) and 
that is partially converted into mechanical work W to lift 
the load. However, according to the Second Law, some 
of this input energy must be dissipated D into the heat 
bath at T\ . If the system has an intrinsic potential there 
will also be a change in the potential energy denoted by 
U. So the law of the total energy conservation is, 



R = D + W + U. 



(9) 



Figure 8: Efficiency r\ as a function of the torque r. 
numerics are done at the optimal regime of the SBM. 



The 



Our interest here is to find the efficiency of the motor, 
which is defined, in general, as ratio between the work 
extracted versus the energy input, 



V 



(R)- 



(10) 



The work W is quite easy to calculate since it is just 
the torque times the angle displaced. Then the work 
performed during a time interval tt — ti can be obtained 



(w) = r(e 1 (t f )-d 1 (t i ))=T(v)(t f -t i ). 



(11) 



The evaluation of (R) needs a much more careful anal- 
ysis. Following [jj, we can write, 



R = 



dV(0 1 (t),0 2 (t)) 

[ 002 



]d0 2 {t). 



(12) 



So by using © , we find that the energy that the bath at 
T 2 transfers to the mechanism in the bath at Ti is, 



R = 



-k(9 1 (t)+^>-e 2 (t))+Vl i (-6 2 (t))}d6 2 (t). (13) 



We have to evaluate this integral for different realizations 
and perform usual statistical averages. We have observed 
that this quantity is always large thus implying impor- 
tant looses of heat. 

To progress further let us first note that the potential 
energy U does not contribute because we have a cyclic 
engine. In fact V 2 {— &2(t)) is periodic and bounded and 
thus it can be discarded in the long time interval studied. 

Thus the final expression for the efficiency that has to 
be evaluated numerically is, 



V 



-(9 1 (t f )~9 l {t l )) 



-mi f (Ut) + <f>-6 2 (t))d0 2 (t)) 



(14) 



In Fig. |H1 the efficiency 77 as a function of the torque r 
is plotted. The specific values we use are around those 
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that were shown to be the optimal ones; Vo = 3, T 2 — 2, 
k = 100, <fi = 0.4 and d = 12. The plot shows a parabolic 
curve for the efficiency with a maximum at the middle 
value of the stall torque. Notice that the maximum of 
the efficiency is still extremely low. Let's also underline 
that at the stall force (r ~ 1-25), the efficiency is zero, 
which means that unavoidable heat transfer occurs 0. 
Therefore, even when (v) — > there's a non-zero heat flux 
transfered from the hotter bath to the cooler one. This 
is the reason why these motors have so low efficiencies. 



III. ANALYTICAL STUDY OF THE CURRENT 
INVERSION 

As we have seen in the numerical results the <f> pa- 
rameter controls a spectacular phenomena. The SBM 
exhibits velocity inversion when this parameter is varied. 
We know that this phenomenon exists and has been stud- 
ied in the literature [l3 |. Our purpose now is to predict 
analytically that this happens in our model. Any kind 
of exact calculation in this system seems impossible and, 
so, some approximations have to be assumed. 

Consider the equations that define our Brownian motor 
@ written in the form, 



f(0 1 )-k(6 1 + 4>-0 2 ) + Si(t), (15) 



0a = -fih) + k{e x + <\> - e 2 ) + &(*), (i6) 

where f(9) is the force exerted by the ratchet potential 
and there is no torque. Let us define now the change of 
variables 113 . 



(17) 



y 



(18) 



The relevant variable x describes the evolution of the cen- 
ter of mass and the "irrelevant" variable y describes the 
relative motion of the two-particle system. Introducing 
a redefinition of the noises, 



£i (*)+&(*) 



&(*)-&(*) 



(19) 



(20) 



and using the fact that y is very small, we can make a 
Taylor expansion up to first order, getting the pair of 
equations, 



x = F(x)+yG(x)+r ll (t), 



y = Q(x) + y(R(x) - 2k) + r, 2 (t). 



(21) 



(22) 
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Figure 9: Comparison of v vs (j> obtained from numerical sim- 
ulations of the model © (circles) and from the analytical ex- 
pression 1231 (solid line). The parameters chosen are d = 4, 
k = 100, Vo = 3 and T 2 = 2T X . 



The explicit expressions for these quantities and further 
mathematical details are in the Appendix. Summarizing, 
we have found that the mean angular velocity v can be 



nguiar velo< 

expressed in the standard form 0> 0, ^| , 

27T , a . 

v<j, = —N(l-eP<f), 



(23) 



where A is a kind of normalization constant that does 
not depend very much on the parameter <fi. The relevant 
quantity /?, which is a function of </>, is given by, 



dx- 



-F(x) 
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Ti 
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This expression can be numerically evaluated using the 
expressions for the different functions that appear in it. 
The value of the parameter N has been estimated by 
using only one numerical value of the velocity obtained 
by the simulation (the minimum value at <p = 0.25). If 
P changes its sign as a function of <f> then the inversion 
phenomenon is predicted. 

In Fig|2|we compare the predicted theoretical depen- 
dence of the velocity as a function of the phase to the 
values obtained by numerical simulations of the model. 
The analytical result, considering the amount of approx- 
imations involved, traces very accurately the current in- 
version phenomenon. The fit is very encouraging and 
enlighting. The positions of the maxima are quite well 
determined, the inversion of current is clearly coincident 
with the simulations and even the little shift near <j) — 
and 4> — ^j- is faithfully reproduced. This means that 
in our approximations we have kept the most dominant 
ingredients. Note that the x-axis goes from to 1.57 
because we have imposed d — 4 (-J = 1.57) in order 
to safely use the analytical expressions derived in the 
Appendix. However, for d = 12 the computation of (3 
predicts qualitatively the current inversion too. 
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Now we will give an intuitive explanation of why the 
motor runs either forward or backwards, depending on 
the value of the parameter <p. The original problem can 
be mapped into a multiplicative noise scenario (see the 
App endix for details) well documented in the literature 
0131110. Then we can straight forward identify in 
our case the effective potential V e ff(x) l|A14jl . which does 
not need to be asymmetric anymore, and the effective 
modulation g e ff(x) (| A13|) . which is responsible for the 
breaking of the symmetries that allow a non-zero mean 
velocity in the motor. What is more, by plotting them 
we can very intuitively see whether and why there is a 
forward or backward flux [Til IT7L ITsf . 

In Fig^Jwe plot, for three significative phase shifts, 
the noise coupling function g e f f (2) and the effective po- 
tential V e ff(x). A careful look at each of them easily 
reveals the current sign in a very intuitive way. The key 
point is the fact that when g e ff(x) is large, the noise ef- 
fects are amplified on one side of the hill of the periodic 
potential and, then, the net current will go in that di- 
rection. Then, it is much clearer that depending on the 
position of g e ff{x) with respect to V e ff(x) the current 
will go forwards, backwards or, in the symmetric case, it 
will cancel. Thus, this simple scheme explains the quali- 
tative behavior of the current inversion. 

Two more comments are to be made. Just for the sake 
of consistency, we have checked two basic conditions that 
must always hold on these kind of motors. The first one 
concerns temperatures; if we set T± — T2, the integral 
vanishes because the denominator becomes constant and, 
therefore, the average velocity is zero. The second test 
is also related with the first one. When we make k tend 
to infinity (this implies a rigid coupling and then only 
one effective temperature), the dependence on x of the 
multiplicative noise disappears again (only additive noise 
is present) and, then, the motor cannot move anywhere 
on average. 



IV. COMMENTS AND CONCLUSIONS 

We have presented and studied a symmetric Brown- 
ian motor. Its efficiency has been numerically obtained 
and other specific properties have been studied. This 
motor has a relevant external parameter (a phase) which 
induces the phenomenon of current inversion. We have 
performed an analytical calculation with appropriate ap- 
proximations to get an expression for the mean velocity 
in terms of the relevant parameters of the model. This 
formal prediction fits very well the data from the numer- 
ical simulations of SBM model. 

An important conclusion is that, albeit the symmet- 
ric motor has a larger efficiency than other mechanical 
Brownian motors, the efficiency of such devices is very 
small. Regardless of the particular properties of these 
kind of heat engines, they are anyhow unrealistic models 
for molecular motors since it is known that these bio- 
logical systems do transform chemical energy into work, 
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Figure 10: Plots of the shape of the effective multiplicative 
noise g E ff(x) 1 A 1311 (dashed line) and the effective potential 
V e ff(x) l|A14|l (solid line) in arbitrary scales. The phase shifts 
are (from top to bottom) cf> — 0.25, (j> = 0.785 and (j> = 1.375, 
which correspond to the fastest negative, zero and fastest pos- 
itive velocities respectively for d — 4. Arrows indicate the 
velocity direction. 



without the intermediate state of burning fuel. Conse- 
quently, one cannot think of these models as realistic ones 
for biological molecular motors. 

Moreover, the mechanical coupling mechanism be- 
tween both baths acts as a very good heat conductor even 
in situations of very small mean velocity. Therefore, the 
efficiency is only a small fraction of that of Carnot. This 
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is in fact a general feature of heat Brownian motors due 
to the fact that in order to rectify thermal fluctuations 
these systems must be put in contact and a lot of heat 
is interchanged. The diode rectifier of Ref. does not 
present this problem and accordingly it can arrive near 
the Carnot efficiency in some limits. 
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Appendix A: ANALYTICAL APPROACH 



92 ^ = 2 I ~2k 



(A9) 



Let's write now the Fokker-Planck equation associated 
to Eq. (|A6|) . 



d t P(x,t) = -d x J(x,t). 



(A10) 



where, 



J{x,t) = H{x)P(x,t) - K B T l \g 1 {x)d x g l {x)P{x,t)] 

-K B T 2 [g 2 (x)d x g 2 (x)P(x,t)}. (All) 

After some manipulations with partial derivatives, the 
probability current above can be rewritten as, 

J(x,t) = H(x)P(x,t) - [g ef f(x)d x g ef f(x)P(x,t)], 

(A12) 

where, 



Our starting point is the pair of coupled equations i|21|) 
and H22f) where the different functions are, 
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(Al) 



(A2) 



(A3) 



(A4) 



Since the variable y has a faster dynamics we will first 
eliminate adiabatically equation l|22|) . and we will also 
discard the function R(x) because the parameter k is 
very large. Then equation (|22|l reduces to, 



(A5) 



Substituting now this expression in equation 1|21|) we get 
a Langevin equation with two multiplicative noises, 



x = H{x) +gi(x)£i(t) + <fe(a06(t), 
where the new functions are, 

H(x) = F(x) + ±Q( x )G(x), 



r ^ 1 (, G 0) 



(A6) 



(A7) 



(A8) 



9eff(x) = \J K B T\g\{x) + K B T 2 g\(x), (A13) 
and H(x) can be related to an effective potential, 

V eff (x) = - [ H(x')dx'. (A14) 



Then we have to solve the equation (|A12|) in the steady 
state for a constant flux J with periodic conditions. The 
first step is to reduce this equation to a Bernoulli form 
which can be formally integrated. By imposing periodic 
boundary conditions, P(x) = P(x + L) (where L = ^), 
we get, 



L e -/3(z) 

J I dx — , 

/O 9eff{x) 



(A15) 



where Pq is a constant that can be found by using the 
normalization condition P{x)dx — 1, and (3 is a rele- 
vant function whose expression is, 



(3{x) 



dx' 



-H(x>) 



fi(x') + (T 2 /T 1 )g*(x'y 



(A16) 



where we take K B T\ = 1. Then the mean velocity v = 
(x) ~ (6*i ) (since (6*i) ~ (9 2 )) is found to be proportional 
to (1 - e^)). What is left to do is to find the (3 integral 
for every 0. 

To simplify the calculation of the integral, we now 
make an expansion in powers of i. Since the value of 
k that makes the motor run faster is around k = 100, 
one can safely suppose that the terms of the order (i) 2 
and so on will not notably contribute to the integral for 
d < 5. Then, the integral (jAlfifl turns out to be J2U, 
which is evaluated numerically. 
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